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Abstract 

We analyze situations where a saddle-node bifurcation occurs on a fractal basin boundary. 
Specifically, we are interested in what happens when a system parameter is slowly swept in 
time through the bifurcation. Such situations are known to be indeterminate in the sense 
that it is difficult to predict the eventual fate of an orbit that tracks the pre-bifurcation node 
attractor as the system parameter is swept through the bifurcation. In this paper we investigate 
the scaling of (1) the fractal basin boundary of the static (i.e., unswept) system near the 
saddle-node bifurcation, (2) the dependence of the orbit's final destination on the sweeping 
rate, (3) the dependence of the time it takes for an attractor to capture a swept orbit on 
the sweeping rate, and (4) the dependence of the final attractor capture probability on the 
noise level. With respect to noise, our main result is that the effect of noise scales with the 
5/6 power of the parameter drift rate. Our approach is to first investigate all these issues 
using one-dimensional map models. The simplification of treatment inherent in one dimension 
greatly facilitates analysis and numerical experiment, aiding us in obtaining the new results 
listed above. Following our one-dimensional investigations, we explain that these results can be 
applied to two-dimensional systems. We show, through numerical experiments on a periodically 
forced second order differential equation example, that the scalings we have found also apply 
to systems that result in two dimensional maps. 
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I. INTRODUCTION 



It is common for dynamical systems to have two or more coexisting attractors. In 
predicting the long-term behavior of a such a system, it is important to determine sets 
of initial conditions of orbits that approach each attractor (i.e., the basins of attrac- 
tion). The boundaries of such sets are often fractal (jl|, Chapter 5 of 0, and references 
therein). The fine-scale fractal structure of such a boundary implies increased sensitivity 
to errors in the initial conditions: Even a considerable decrease in the uncertainty of 
initial conditions may yield only a relatively small decrease in the probability of mak- 
ing an error in determining in which basin such an initial condition belongs 0, 0| • For 
discussion of fractal basin boundaries in experiments, see Chapter 14 of 0. 

Thompson and Soliman 0] showed that another source of uncertainty induced by 
fractal basin boundaries may arise in situations in which there is slow (adiabatic) varia- 
tion of the system. For example, consider a fixed point attractor of a map (a node). As 
a system parameter varies slowly, an orbit initially placed on the node attractor moves 
with time, closely following the location of the solution for the fixed point in the absence 
of the temporal parameter variation. As the parameter varies, the node attractor may 
suffer a saddle-node bifurcation. For definiteness, say that the node attractor exists for 
values of the parameter /i in the range /i < /x*, and that the saddle- node bifurcation of 
the node occurs at n — /i*. Now assume that, for a parameter interval with 
fiL < /i* < fiR, in addition to the node, there are also two other attractors A and B, 
and that the boundary of the basin of attractor A, attractor B and the node is a fractal 
basin boundary. We are interested in the typical case where, before the bifurcation, 
the saddle lies on the fractal basin boundary, and thus, at the bifurcation, the merged 
saddle-node orbit is on the basin boundary. In such a case an arbitrarily small ball about 
the saddle-node at /i = /i* contains pieces of the basins of both A and B. Thus, as ji 
slowly increases through /i*, it is unclear whether the orbit following the node will go to 
A or to B after the node attractor is destroyed by the bifurcation. In practice, noise or 
round-off error may lead the orbit to go to one attractor or the other, and the result can 
often depend very sensitively on the specific value of the slow rate at which the system 
parameter varies. 

We note that the study of orbits swept through an indeterminate saddle-node bifurca- 
tion belongs to the theory of dynamical bifurcations. Many authors have analyzed orbits 
swept through other bifurcations, like the period doubling bifurcation p, the pitchfork 
bifurcation p, 0], and the transcritical bifurcation Q. In all these studies of the bifur- 
cations listed above, the local structure before and after the bifurcation includes stable 
invariant manifolds varying smoothly with the bifurcation parameter (i.e., a stable fixed 
point that exists before or after the bifurcation, and whose location varies smoothly with 
the bifurcation parameter). This particular feature of the local bifurcation structure, not 
shared by the saddle-node bifurcation, allows for well-posed, locally defined, problems 
of dynamical bifurcations. The static saddle-node bifurcation has received much atten- 
tion in theory and experiments (USE); but so far, no dynamical bifurcation problems 
have been defined for the saddle-node bifurcation. In this work, we demonstrate that, 
in certain common situations, global structure (i.e., an invariant Cantor set or a fractal 
basin boundary) adds to the local properties of the saddle-node bifurcation and allows 
for well-posed problems of dynamical bifurcations. 

Situations where a saddle-node bifurcation occurs on a fractal basin boundary have 
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been studied in two dimensional Poincare maps of damped forced oscillators 0, E, 12l . 
Several examples of such systems are known 0, 2|> and it seems that this is a common 
occurence in dynamical systems. In this work, we first focus on saddle-node bifurcations 
that occur for one parameter families of smooth one dimensional maps having multiple 
critical points (a critical point is a point at which the derivative of the map vanishes). 
Since one dimensional dynamics is simpler than two dimensional dynamics, indeterminate 
bifurcations can be more simply studied, without the distraction of extra mathematical 
structure. Taking advantage of this, we are able to efficiently investigate several scaling 
properties of these bifurcations. In particular, we investigate the scaling of (1) the fractal 
basin boundary of the static (i.e., unswept) system near the saddle-node bifurcation 
( Sees. ITTBl and lrTT^) . (2) the dependence of the orbit final destination on the sweeping rate 
(Sec. Ill D|) . (3) the dependence of the time it takes for an attractor to capture the swept 
orbit following the bifurcation on the sweeping rate fSec. I11FJ|) . and (4) the dependence 
of the final attractor capture probability on the noise level (Sec. IllFj) . Following our 
one-dimensional investigations, we explain that these results apply to two dimensional 
systems. We show, through numerical experiments on the periodically forced Duffing 
oscillator, that the scalings we have found also apply to higher dimensional systems 
(Sec. EH. 

For one-dimensional maps, a situation dynamically similar to that in which there is 
indeterminacy in which attractor captures the orbit can also occur in cases where there 
are two rather than three (or more) attractors (Sec. I1VJ) . In particular, we can have 
the situation where one attractor persists for all values of the parameters we consider, 
and the other attractor is a node which is destroyed via a saddle-node bifurcation on 
the basin boundary separating the basins of the two attractors. In such a situation, an 
orbit starting on the node, and swept through the saddle-node bifurcation, will go to the 
remaining attractor. It is possible to distinguish different ways that the orbit initially 
on the node approaches the remaining attractor. We find that the way in which this 
attractor is approached can be indeterminate. 



II. INDETERMINACY IN WHICH ATTRACTOR IS APPROACHED 

We consider the general situation of a one dimensional real map fn(x) depending on 
a parameter /x. We assume the following: (1) the map is twice differentiable with respect 
to x, and once differentiable with respect to /x (the derivatives are continuous); (2) f ^ 
has at least two attractors sharing a fractal basin boundary for parameter values in the 
vicinity of /x*; and (3) an attracting fixed point 2* of the map f^x) is destroyed by a 
saddle-node bifurcation as the parameter /x increases through a critical value /x*, and 
this saddle-node bifurcation occurs on the common boundary of the basins of the two 
attractors. 

We first recall the saddle-node bifurcation theorem (see for example [H). If the map 
fn(x) satisfies: (a) /^(x*) = x*, (b) ^-(ar*) = 1, (c) ^fi^{x*) > 0, and (d) §£(x*;/x*) > 
0, then the map undergoes a backward saddle- node bifurcation (i.e., the node attractor 
is destroyed at x* as /x increases through /x*). If the inequality in either (c) or (d) is 
reversed, then the map undergoes a forward saddle-node bifurcation, while, if both these 
inequalities are reversed, the bifurcation remains backward. A saddle-node bifurcation 
in a one dimensional map is also called a tangent or a fold bifurcation. 
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A. Model 



As an illustration of an indeterminate saddle-node bifurcation in a one-dimensional 
map, we construct an example in the following way. We consider the logistic map for 
a parameter value where there is a stable period three orbit. We denote this map g(x) 
and its third iterate g^(x). The map g^\x) has three stable fixed points. We perturb 
the map g^(x) by adding a function (which depends on a parameter /i) that will cause 
a saddle-node bifurcation of one of the attracting fixed points but not of the other two 
[see Figs. 1(a) and 1(b)]. We investigate 

fn{x) = 9^{x) + / usin(37rx), where g(x) = 3.832 x(l — x). (1) 

Numerical calculations show that the function f^(x) satisfies all the conditions of the 
saddle-node bifurcation theorem for having a backward saddle-node bifurcation at ~ 
0.15970 and /i* ~ 0.00279. Figure 2(a) displays how the basins of the three attracting 
fixed points of the map f ^ change with variation of fi. For /i = the third iterate of 
the logistic map is unperturbed, and it has three attracting fixed points whose basins we 
color-coded with blue, green and red. For every value of /i, the red region R\p] is the set 
of initial conditions attracted to the rightmost stable fixed point which we denote R^. 
The green region G[/j] is the set of initial conditions attracted to the middle stable fixed 
point which we denote G M . The blue region B\p] is the set of initial conditions attracted 
to the leftmost stable fixed point which we denote _B M . 

For each of these colored sets has infinitely many disjoint intervals and a 

fractal boundary. As /i increases, the leftmost stable fixed point is destroyed via a 
saddle-node bifurcation on the fractal basin boundary. In fact, in this case, for fi < /i*, 
every boundary point of one basin is a boundary point for all three basins. (That is, an 
arbitrarily small x-interval centered about any point on the boundary of any one of the 
basins contains pieces of the other two basins.) The basins are so-called Wada basins 
This phenomenon of a saddle-node bifurcation on the fractal boundary of Wada basins 
also occurs for the damped forced oscillators studied in Refs. 11, Alternatively, if we 



look at the saddle- node bifurcation as \i decreases through the value //*, then the basin 
B\p] of the newly created stable fixed point immediately has infinitely many disjoint 
intervals and its boundary displays fractal structure. According to the terminology of 



Robert et al. |lj|, we may consider this bifurcation an example of an 'explosion'. 



B. Dimension of the Fractal Basin Boundary 

Figure 3 graphs the computed dimension D of the fractal basin boundary versus the 
parameter \i. For /i < we observe that D appears to be a continuous function of 
fi. Park et al. ^| argue that the fractal dimension of the basin boundary near /x*, for 
/j < /i*, scales as 

D(n) « D, - k{fj„ - /x) 1/2 , (2) 

with the dimension at fi = //* (D* is less than the dimension of the phase space), 
and k a positive constant. Figure 3 shows that the boundary dimension D experiences 
a discontinuous jump at the saddle-node bifurcation when \x = /i*. We believe that this 
is due to the fact that the basin B[/j] suddenly disappears for /i > 

The existence of a fractal basin boundary has important practical consequences. In 
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particular, for the purpose of determining which attractor eventually captures a given 
orbit, the arbitrarily fine-scaled structure of fractal basin boundaries implies considerable 
sensitivity to small errors in initial conditions. If we assume that initial points cannot be 
located more precisely than some e > 0, then we cannot determine which basin a point 
is in, if it is within e of the basin boundary. Such points are called e-uncertain. The 
Lebesgue measure of the set of e-uncertain points (in a bounded region of interest) scales 
like e D °~ D , where D is the dimension of the phase space (D = 1 for one dimensional 
maps) and D is the box-counting dimension of the basin boundary For the case of a 
fractal basin boundary (Dq — D) < 1. When Dq — D is small, a large decrease in e results 
in a relatively small decrease in e Do ~ D . This is discussed in Ref. ^ which defines the 
uncertainty dimension, D u , as follows. Say we randomly pick an initial condition x with 
uniform probability density in a state-space region S. Then we randomly pick another 
initial condition y in S, such that \y — x\ < e. Let p(e, S) be the probability that x and 
y are in different basins. [We can think of p(e, S) as the probability that an error will be 
made in determing the basin of an initial condition if the initial condition has uncertainty 
of size e.] The uncertainty dimension of the basin boundary D u is defined as the limit 
of lnp(e, S)/ ln(e) as e goes to zero Thus, the probability of error scales as p(e, S) ~ 
e D °~ D ", where for fractal basin boundaries Dq — D u < 1. This indicates enhanced 
sensitivity to small uncertainty in initial conditions. For example, if Dq — D u = 0.2, then 
a decrease of the initial condition uncertainty e by a factor of 10 leads to only a relative 
small decrease in the final state uncertainty p(e, S), since p decreases by a factor of about 
10 2 ~ 1.6. Thus, in practical terms, it may be essentially impossible to significantly 
reduce the final state uncertainty. In Ref. |l| it was conjectured that the box-counting 
dimension equals the uncertainty dimension for basin boundaries in typical dynamical 
systems. In Ref. it is proven that the box-counting dimension, the uncertainty 
dimension and the Hausdorff dimension are all equal for the basin boundaries of one and 
two dimensional systems that are uniformly hyperbolic on their basin boundary. 

We now explain some aspects of the character of the dependence of D on /i (see Fig. 3). 
From Refs. [l8| it follows that the box-counting dimension and the Hausdorff dimension 
coincide for all intervals of /i for which the map f ^ is hyperbolic on the basin boundary, 
and that the dimension depends continuously on the parameter [i in these intervals. 
For there are many parameter values for which the map has a saddle- node 

bifurcation of a periodic orbit on the fractal basin boundary. At such parameter values, 
which we refer to as saddle-node bifurcation parameter values, the dimension is expected 
to be discontinuous (as it is at the saddle- node bifurcation of the fixed point, = 
see Fig. 3). In fact, there exist sequences of saddle-node bifurcation parameter values 
converging to /z* [l6j. Furthermore, for each parameter value /i > /i* for which the map 
undergoes a saddle-node bifurcation, there exists a sequence of saddle-node bifurcation 
parameter values converging to that parameter value. The basins of attraction of the 
periodic orbits created by saddle-node bifurcations of high period exist only for very 
small intervals of the parameter /i. We did not encounter them numerically by iterating 
initial conditions for a discrete set of values of the parameter /z, as we did for the basin 
of our fixed point attractor. 
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C. Scaling of the Fractal Basin Boundary 



Just past /i*, the remaining green and red basins display an alternating stripe struc- 
ture [see Fig. 2(b)]. The red and green stripes are interlaced in a fractal structure. As we 
approach the bifurcation point, the interlacing becomes finer and finer scaled, with the 
scale approaching zero as fi approaches /i*. Similar fine scaled structure is present in the 
neighborhood of all preiterates of x„. If one changes the horizontal axis of Figs. 2(a,b) 
from fi to (fi — ii*)" 1 / 2 , then, the complex alternating stripe structure appears asymp- 
totically periodic [see Fig. 4(a)]. [Thus, with identical horizontal scale, the dimension 
plot in Fig. 4(b) appears asymptotically periodic, as well.] We now explain why this is 
so. We restrict our discussion to a small neighborhood of x*. Consider the second order 
expansion of in the vicinity of x* and //* 



fa(x) = u + x + ax 2 , where l X X X *' (3) 
^ v ' \n = fi-fi*, v ' 

and a ~ 89.4315. The trajectories of fp, in the neighborhood of x — 0, for jl close to zero, 
are good approximations to trajectories of f^ in the neighborhood of x = x*, for fi close 
to /i*. Assume that we start with a certain initial condition for fp, and we ask 

the following question: What are all the positive values of the parameter fi such that a 
trajectory passes through a fixed position Xf > at some iterate n? For any given Xf 
which is not on the fractal basin boundary, there exists a range of ft such that iterates 
of Xf under f^ evolve to the same final attractor, for all values of fi in that range. In 
particular, once ax 2 appreciably exceeds fi, the subse quent evolution is approximately 
independent of fi. Thus, we can choose x/ ^> a/ fi/a, but still small enough so that it 
lies in the region of validity of the canonical form ©• There exists a range of such Xf 
values satisfying these requirements provided that \fl\ is small enough. 

Since consecutive iterates of fp, in the neighborhood of x = for fi close to zero differ 
only slightly, we approximate the one dimensional map, 

x n+1 = fp(x n ) = fi + x n + axl, (4) 

by the differential equation 0, 

— = fi + ax , (5) 
an 

where in (j3J) n is considered as a continuous, rather than a discrete, variable. Integrating 
(jSJ) from x s to £f yields 

n^Jafi = arctan ^J^Xfj — arctan (^-^^^^j ■ (6) 

Close to the saddle-node bifurcation (i.e., < /I < 1, and x s j close to zero), fp is a 
good approximation to / M . For \x s j\y/ (a/fi) ^> 1 Eq. © becomes 

n\fafi ~ 7T. (7) 

jt^ rest 

The values of fi n satisfying Eq. (JJj) increase with n in step of \fa/n. For our example 
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we have a ~ 89.4315, thus \fa/ir pa 3.010. Counting many periods like those in Fig. 4 in 
the region of x c , the closest critical point to x„ [see Fig. 4(a)], we find that the period of 
the stripe structure is 3.015, which is in good agreement with our theoretical value. 

In order to investigate the structure of the fractal basin boundary in the vicinity of 
the saddle-node bifurcation (i.e., x s close to £* = 0), we consider © in the case where 
we demand only \xf\\J{a/ p) ^> 1. Thus, Eq. (0) becomes 

71 ( 
n v au ~ arctan , / — x s 

2 VV V 

Let fin (x s ) denote the solution of Eq. (JHJ) for ft,. Equation (jHJ) implies the behavior 
of fin^^ixs) as function of x s and n as sketched in Fig. 5. For a fixed n, fin 1 ^ 2 has 
a horizontal asymptote at the value ny/a/n as x s — > — oo, and a vertical asymptote to 
infinity at x s = l/(an). For x s < 0, we have an infinite number of values of the parameter 
fi, for which an orbit of starting at x s passes through the same position Xf, after some 
number of iterations. For x s = (i.e., x s = x*), we also have an infinite number of 
fi n (0), but with constant step 2y/a/n rather than y/a/n (see the intersections marked 



1/2 



with black dots in the Fig. 5). This is hard to verify from numerics, since — (0) = 

a 3 ^ 2 (2n/ir) 2 increases with n 2 , and the stripes become very tilted in the neighborhood 
of x s = x* = 0. [See Fig. 4(a), where the approximate positions of x c and x* on the 
vertical axis are indicated.] For x s > 0, \x n has only a limited number of values with 

^max < 1/(OX ). 



D. Sweeping Through an Indeterminate Saddle-Node Bifurcation 

In order to understand the consequences of a saddle-node bifurcation on a fractal basin 
boundary for systems experiencing slow drift, we imagine the following experiment. We 
start with the dynamical system at parameter /i s < //*, with x on the attractor to 
be destroyed at fi = //* by a saddle- node bifurcation (i.e., B u ). Then, as we iterate, we 
slowly change /i by a small constant amount 5/x per iterate, thus increasing \x from /i s to 
fjif > /i*, 

[In = yUs +nS/ji. 

When /i > fif we stop sweeping the parameter /i, and, by iterating further, we determine 
to which of the remaining attractors of f Uf the orbit goes. Numerically, we observe 
that, if (jif — Ai*) is not too small, then, by the time ///is reached, the orbit is close 
to the attractor of f Uf to which it goes. [From our subsequent analysis, 'not too small 
|a*s,/ ~ A** I' translates to choices of 5/j, that satisfy (Sji) 2 ^ 3 \fx s j — We repeat this 
for different values of S/j, and we graph the final attractor position for the orbit versus 
5/i [see Fig. 6(a)]. For convenience in the graphical representation of Figs. 6(a,b), we 
have represented the attractor of the green region G\p], denoted G Uf , as a 0, and the 
attractor of the red region R\p], denoted R Uf , as a 1. In Fig. 6(a) we use of 25,000 points 
having the vertical coordinate either or 1, which we connect with straight lines. In 
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an interval of 8fi for which the system reaches the same final attractor (either or 1), 
the lines connecting the points are horizontal. Such intervals appear as white bands in 
Fig. 6, if they are wider than the width of the plotted lines connecting O's and l's. For 
example, in Fig 6(a), the white band centered at 5fi = 0.8 x 10~ 3 has at the bottom a 
thick horizontal line, which indicates that for the whole of that interval, the orbit reaches 
the attractor G fJiJ which we represented by 0. Adjacent intervals of width less than the 
plotted lines appear as black bands. Within such black bands, an uncertainty in £/i of 
size equal to the width of the plotted line makes the attractor that the orbit goes to 
indeterminate. Figure 6(a) shows that the widths of the white bands decrease as S/j, 
decreases, such that, for small we see only black. 

If (/if — /i*) is large enough (i.e., (<5/i) 2 / 3 <C — //*|), numerics and our subsequent 
analysis show that Fig. 6 is independent of /i/. This fact can be understood as follows. 
Once fi = fif, the orbit typically lands in the green or the red basin of attraction and goes 
to the corresponding attractor. Due to sweeping, it is possible for the orbit to switch 
from being in one basin of attraction of the time-independent map f ^ to the other, since 
the basin boundary between G\jj] and R\p] changes with fi. However, the sweeping of \i 
is slow (i.e., 5fi is small), and, once (// — //*) is large enough, the orbit is far enough from 
the fractal basin boundary, and the fractal basin boundary changes too little to switch 
the orbit between G[/j] and R[jj\. 

We also find numerically that Figs. 6(a,b) are independent of the initial condition 
Xq, provided that it is in the blue basin B[fi s ], sufficiently far from the fractal basin 
boundary, and that \/j s — is not too small (i.e., (<5/i) 2 / 3 <C \/j s — //*|). 

If one changes the horizontal scale of Fig. 6(a) from <5/i to 1/5 /i [see Fig. 6(b)], the 
complex band structure appears asymptotically periodic. Furthermore, we find that the 
period in (1/5/x) of the structure in Fig. 6(b) asymptotically approaches —l/(fi s — /x*) 
as 8 ii becomes small. 

In order to explain this result, we again consider the map f^, the local approximation 
of ffj, in the region of the saddle-node bifurcation. Equations can be approximated 
by 

x n+l — fp, n {%n) — An + °^n' (10) 

fin = frs + n8/2. 

We perform the following numerical experiment. We consider orbits of our approximate 
two dimensional map given by Eq. (fTUj) starting at x s = —^J—p, s /a. We define a final 
state function of an orbit swept with parameter 5fi in the following way. It is if 
the orbit has at least one iterate in a specified fixed interval far from the saddle-node 
bifurcation, and is 1, otherwise. In particular, we take the final state of a swept orbit 
to be if there exists n such that 100 < x n < 250, and to be 1 otherwise. Figure 6(c) 
graphs the corresponding numerical results. Similar to Fig. 6(b), we observe periodic 
behavior in l/5fi with period —l/fi s . In contrast to Fig. 6(b) where the white band 
structure seems fractal, the structure within each period in Fig. 6(c) consists of only 
one interval where the final state is and one interval where the final state is 1. This is 
because 100 < x < 250 is a single interval, while the green basin [denoted in Fig. 6(b)] 
has an infinite number of disjoint intervals and a fractal boundary (see Fig. 2). 

With the similarity between Figs. 6(b) and 6(c) as a guide, we are now in a position 
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to give a theoretical analysis explaining the observed periodicity in 1/ S/i. In particular, 
we now know that this can be explained using the canonical map (JTOJ), and that the 
periodicity result is thus universal [i.e., independent of the details of our particular 
example, Eq. (JTJ)]. For slow sweeping (i.e., S/j small), consecutive iterates of (fTU|) in the 
vicinity of x = and ft = differ only slightly, and we further approximate the system 
by the following Ricatti differential equation, 

dec 

— = ft s + nSu + ax 2 . (11) 
an 

The solution of Eq. (fTT]) can be expressed in terms of the Airy functions Ai and Bi and 
their derivatives, denoted by Ai' and Bi f , 



V Ai'(0 + Bt'(0 (S/j 

,2 



1/3 



f(n)= vm+m • (12) 



where 



1/3 P-s + nS/j 



«"' = -°'l^r- < 13 > 

and rj is a constant to be determined from the initial condition. We are only interested 
in the case of slow sweeping, 8/j <C 1, and x(0) = x s = —a/— /is/a (which is the stable 
fixed point of destroyed by the saddle- node bifurcation at ft — 0). In particular, 
we will consider the case where ft s < and \ft s \ ^> 5/i 2 / 3 (i.e., |^(0)| ^> 1). Using 
x(0) = —yJ—fL s /a to solve for rj yields rj ~ O[^(0)e 2 ^} ^> 1. For positive large values 
of £(n) (i.e., for n small enough), using the corresponding asymptotic expansions of the 
Airy functions [l^, the lowest order in 5fi approximation to (fT^j) is 



4(n)s3 _y_^±^, (14 ) 

with the correction term of higher order in S/j, being negative. Thus, for n sufficiently 
smaller than —ft s /S/M, the swept orbit lags closely behind the fixed point for with ft 
constant. For £ < 0, we use the fact that rj is large to approximate (fT^jl as 

Note that 



gives the lag of the swept orbit relative to the fixed point attractor evaluated at the 
saddle-node bifurcation. Equation (fT5j) does not apply for n > n max , where n max is the 
value of n for which £(n max ) = £, the largest root of Ai(£) = (i.e., £ = —2.3381...). 
At n = n max , the normal form approximation predicts that the orbit diverges to +oo. 
Thus, for n near n max , the normal form approximation of the dynamical system ceases 
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to be valid. Note, however, that (JT3J) can be valid even for £(n) close to £(n max ). This is 
possible because Sfi is small. In particular, we can consider times up to the time n' where 
n' is determined by £' = £,(n') — £ + 5£, (5£ > is small,) provided |x(n')| <C 1 so that 
the normal form applies. That is, we require [Ai' / Ai(£')](6 (i/ a 2 ) 1 / 3 <C 1, which can 
be satisfied even if [Ai'(£')/Ai(£')] is large. Furthermore, we will take the small quantity 
5£ to be not too small (i.e., 5^/(a5fx) 1 ^ 3 ^> 1), so that (n ma x _n ') ^ 1- We then consider 
(J15)) in the range, —(fi s /5fi) < n < n', where the normal form is still valid. 

We use Eq. (fT3j) for answering the following question: What are all the values of the 
parameter 5{i (5{i small) for which an orbit passes exactly through the same position Xf > 
0, at some iterate np. All such orbits would further evolve to the same final attractor, 
independent of 5/i, provided ax 2 ^> £i s +rif S/i; i.e., Xf is large enough that fif = fi s +nf 5/i 
does not much influence the orbit after x reaches Xf. [Denote £(%) as = Using 
(fTol we can estimate when this occurs, ax 2 = [Ai'(£f)/ Ai(^f)] 2 (5fx 2 /a) 1 / 3 (jl s + nf 5fj>) 

or [Ai 1 '(£/) / Ai^f)] 2 ^> £/. This inequality is satisfied when gets near £, which is the 
largest zero of Ai (i.e., = £ + 5£, where <5£ is a small positive quantity). We now 
rewrite Eq. (J15)) in the following way 



representing a transcedental equation in 5fi where fi s and Xf are fixed, rif is a large 
positive integer (i.e., n/ — 1 is the integer part of {/if — fi s )/5fi), and K(Q is the in- 
verse function of Ai'(£)/Ai(£) in the neighborhood of £ = (a 2 /5/i) 1//3 xj 1. Thus 



K[{a 2 /5fi) l / 3 x f }\ < \K(oo)\ = |||. The difference [l/5fi{x f ,n f + 1) - l/5/i(x f ,n f )}, 



where 5/i(xf,nf) is the solution of Eq. (fTTj) . yields the limit period of the attracting state 
versus l/Sfi graph (see Fig. 6). We denote this limit period by A (l/<5/z). For small 
the term involving K[(a 2 /Sfi) 1 ^ 3 xf\ in Eq. (jl7j) can be neglected, and we get A (1/5/x) = 
— (i^ 1 = (—fi s + /i*) -1 . Figure 7 graphs numerical results of [A (1/5/i)] -1 versus /i s for 
our map example given by Eq. ©. The fit line is [A (1/5/i)]" 1 = -0.9986/i s + 0.0028, 
which agrees well with the prediction of the above analysis and our numerical value for 
fi at the bifurcation, /i* « 0.00279. 

An alternate point of view on this scaling property is as follows. For /i < (i.e., 
ji < fiA and slow sweeping (i.e., <5/i small), the orbit closely follows the stable fixed 
point attractor of until /x > 0, and the saddle-node bifurcation takes place. However, 
due to the discreteness of n, the first nonnegative value of /i depends on fi s and 5li (see 
Fig. 8). Now consider two values of 5li, one 5ii m satisfying fi s + m5fi m = 0, and another 
<5/i m+ i satisfying jl s + (m + 1) <5/i m+ i = 0. Because 5\i m and £/i m+ i are very close (for 
large m) and both lead /i{n) to pass through /t = /2* = (one at time n — m, and the 
other at time n — m+ 1), it is reasonable to assume that their orbits for /l s /5ll < n <n' 
are similar (except for a time shift n — > n + 1); i.e., they go to the same attractor. Thus, 
the period of is approximately A (1/5ll) = l/5fi m+ i — l/5ii m = — /C 1 - 

We now consider the intervals of 1/5 fx between the centers of consecutive wide white 
bands in Fig. 6(b). Figure 9 graphs the calculated fractal dimension D' of the boundary 
between white bands in these consecutive intervals versus their center value of 1/5 fx. 
From Fig. 9, we see that as l/5fi increases, the graph of the fractal dimension D' does 
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(17) 
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not converge to a definite value, but displays further structure. Nevertheless, numerics 
show that as l/Sfi becomes large (i.e., in the range of 6.5 x 10 5 ), D' varies around the 
value 0.952. This is consistent with the numerics presented in Fig. 4(b) which graphs 
the dimension of the fractal basin boundary for the time-independent map / M , at fixed 
values of the parameter // where //>//*. Thus, for large D' provides an esti- 

mate of the dimension of the fractal basin boundary in the absence of sweeping at /i > /i*. 

We now discuss a possible experimental application of our analysis. The conceptually 
most straightforward method of measuring a fractal basin boundary would be to repeat 
many experiments each with precisely chosen initial conditions. By determining the final 
attractor corresponding to each initial condition, basins of attraction could conceivably 
be mapped out |3|. However, it is commonly the case that accurate control of initial 
conditions is not feasible for experiments. Thus, the application of this direct method is 
limited, and, as a consequence, fractal basin boundaries have received little experimental 
study, in spite of their fundamental importance. If a saddle-node bifurcation occurs on 
the fractal basin boundary, an experiment can be arranged to take advantage of this. 
In this case, the purpose of the experiment would be to measure the dimension D' 
as an estimate of the fractal dimension of the basin boundary D. The measurements 
would determine the final attractor of orbits starting at the attractor to be destroyed by 
the saddle-node bifurcation, and swept through the saddle-node bifurcation at different 
velocities (i.e., the experimental data corresponding to the numerics in Fig. 6). This does 
not require precise control of the initial conditions of the orbits. It is sufficient for the 
initial condition to be in the basin of the attractor to be destroyed by the saddle-node 
bifurcation; after enough time, the orbit will be as close to the attractor as the noise 
level allows. Then, the orbit may be swept through the saddle-node bifurcation. The 
final states of the orbits are attractors; in their final states, orbits are robust to noise and 
to measurement perturbations. The only parameters which require rigorous control are 
the sweeping velocity (i.e., and the initial value of the parameter to be swept (i.e., 
/i s ); precise knowledge of the parameter value where the saddle-node bifurcation takes 
place (i.e., is not needed. [It is also required that the noise level be sufficiently low 
(see Sec. EE)-] 



E. Capture Time 



A question of interest is how much time it takes for a swept orbit to reach the final 
attracting state. Namely, we ask how many iterations with \i > /i* are needed for the 
orbit to reach a neighborhood of the attractor having the green basin. Due to slow 
sweeping, the location of the attractor changes slightly on every iterate. If x^ is a fixed 
point attractor of / M (with fi constant), then a small change 5fi in the parameter /i, yields 
a change in the position of the fixed point attractor, 



\%H+5fj, ~ x n) = Sx — 5fl qj^ 



We consider the swept orbit to have reached its final attractor if consecutive iterates differ 
by about Sx (which is proportional to 5fi). For numerical purposes, we consider that the 
orbit has reached its final state if \x n+ i — x n \ < 105//. In our numerical experiments, 
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this condition is satisfied by every orbit before \i reaches its final value ///. We refer 
to the number of iterations with needed to reach the final state as the capture 

time of the corresponding orbit. Figure 10 plots the capture time by the attractor G M/ 
[having the green basin in Fig. 2] versus 1/5/i for a range corresponding to one period 
of the structure in Fig. 6(b). No points are plotted for values of 6fi for which the orbit 
reaches the attractor R flf . The capture time graph has fractal features, since for many 
values of 5 Li the orbit gets close to the fractal boundary between R\jj] and G[fi}. Using 
the fact that the final destination of the orbit versus 1/Sli is asymptotically periodic [see 
Fig. 6(b)], we can provide a further description of the capture time graph. We consider 
the series of the largest intervals of 1/5 fi for which the orbit reaches the attractor G llj . 
[see Fig. 6(b); we refer to the wide white band around 1/Sfi = 2400 and the similar ones 
which are (asymptotically) separated by an integer number of periods]. Orbits swept 
with 5fi at the centers of these intervals spend only a small number of iterations close to 
the common fractal boundary of R\p\ and G[fi]. Thus, the capture time of such similar 
orbits does not depend on the structure of the fractal basin boundary. We use Eq. (fTo} 
as an approximate description of these orbits. A swept orbit reaches its final attracting 
state as x{n) becomes large. Then, the orbit is rapidly trapped in the neighborhood of 
one of the swept attractors of f^. Thus, we equate the argument of the Airy function in 
the denominator to its first root [see (|T5j)]. solve for n, and substract —£i s /5li (the time 
for fx to reach the bifurcation value). This yields the following approximate formula for 
the capture time 

nc^lfl(a^)- 1 / 3 , (18) 

where £ = —2.3381... is the largest root of the Airy function Ai. Thus, we predict that for 
small a log-log plot of the capture time of the selected orbits versus 6fi is a straight 
line with slope -1/3. Figure 11 shows the corresponding numerical results. The best 
fitting line (not shown) has slope -0.31, in agreement with our prediction j20j. 



F. Sweeping Through an Indeterminate Saddle-Node Bifurcation in the Pres- 
ence of Noise 

We now consider the addition of noise. Thus, we change our swept dynamical system 

to 

x n +i = Un( x n) + Ae n , (19) 

where e n is random with uniform probability density in the interval [—1, 1], and A is a 
parameter which we call the noise amplitude. See Fig. 6(a) which shows the numerical 
results of the final destination of the orbits versus 5/i in the case A = 0. The graph ex- 
hibits fractal features of structure at arbitrarily small scales. The addition of small noise 
is expected to alter this structure, switching the final destination of orbits. In this case, 
it is appropriate to study the probability of orbits reaching one of the final destinations. 
For every A, we compute the final attractor of a large number of orbits having identical 
initial condition and parameters, but with different realizations of the noise. We estimate 
the probability that an orbit reaches a certain attractor by the fraction of such orbits 
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that have reached the specified attractor in our numerical simulation. Figure 12 graphs 
the probability that an orbit reaches the attractor G M/ versus the noise amplitude A. We 
present five graphs corresponding to five different values of Sfi equally spaced in a range of 
10" 7 centered at 1(T 5 (i.e., 5^ = 10" 5 , 10~ 5 ±2.5 x 10~ 8 and 10" 5 ±5 x 10~ 8 ). We notice 
that the probability graphs have different shapes, but a common horizontal asymptote 
in the limit of large noise. The value of the horizontal asymptote, approximately equal 
to 0.5, is related to the relative measure of the corresponding basin. 

As in the previous subsection, we take advantage of the asymptotically periodic struc- 
ture of the noiseless final destination graph versus 1/5 fj, [see Fig. 6(b)]. We consider 
centers of the largest intervals of l/5fi for which an orbit reaches the middle attractor in 
the absence of noise. We chose five such values of 5fi, spread over two decades, where the 
ratio of consecutive values is approximately 3. Figure 13(a) graphs the probability that 
an orbit reaches the middle fixed point attractor versus the noise amplitude A, for the 
five selected values of 5fi. We notice that all the curves have qualitatively similar shape. 
For a range from zero to small A, the probability is 1, and as A increases, the probability 
decreases to a horizontal asymptote. The rightmost curve in the family corresponds to 
the largest value of 5fi (<5/i ~ 3.445974 x 10~ 5 ), and the leftmost curve corresponds to 
the smallest value of 5fi (5fi ~ 4.243522 x 10~ 7 ). Figure 13(b) shows the same family of 
curves as in Fig. 13(a), but with the horizontal scale changed from A to A/(5fj,) 5 ^ 6 . All 
data collapse to a single curve, indicating that the probability that a swept orbit reaches 
the attractor G flf depends only on the reduced variable A/(5fi) 5 ^ 6 . Later, we provide a 
theoretical argument for this scaling. 

In order to gain some understanding of this result, we follow the idea of Sec. Ill 131 
and use the canonical form fa to propose a simplified setup of our problem. We 
modify (fTUj) by the addition of a noise term Ae n in the right hand side of the first 
equation of (fTUj) . We are interested in the probability that a swept orbit has at least 
one iterate, specified fixed interval far from the vicinity of the saddle-node 

bifurcation. More precisely, we analyze how this probability changes versus A and 
5fi. Depending on the choice of interval and the choice of 5fi, the probability versus 
A graph (not shown) has various shapes. For numerical purposes, we choose our 
fixed interval to be the same as that of Sec. Ill Ul 100 < x < 250. We then select 
values of 5fi for which a noiseless swept orbit, starting at x s = —y/—ji s /a, reaches 
exactly the center of our fixed interval. The inverse of these values of 5fi are centers 
of intervals where the final state of the swept orbits is [see Fig. 6(c)]. We consider 
five such values of 5fi, where the ratio of consecutive values is approximately 3. Figure 
14(a) shows the probability that a swept orbit has an iterate in our fixed interval 
versus the noise amplitude for the selected values of 5fi. Figure 14(a) shares the 
qualitative characteristics of Fig. 13(a), with the only noticeable difference that the 
value of the horizontal asymptote is now approximately 0.1. Figure 14(b) shows the 
same family of curves as in Fig. 14(a), where the horizontal scale has been changed 
from A to A/(5/j) 5 ^ 6 . As for Fig. 12(b), this achieves good collapse of the family of curves. 

We now present a theoretical argument for why the probability of reaching an attractor 
depends on <5/i and A only through the scaled variable A/(5fx) 5 ^ 6 when 5fi and A are 
small. From our results in Figs. 14, we know that the scaling we wish to demonstrate 
should be obtainable by use of the canonical form fa. Accordingly, we again use the 
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differential equation approximation (jllj) . but with a noise term added, 

doc 

— = n5fi + ax 2 + Ae(n) , (20) 
an 

where e(n) is white noise, 

(e»> = 0, (e(n + n»)> = 

and we have redefined the origin of the time variable n so that the parameter fx sweeps 
through zero at n = (i.e., we replaced n by n— \(1 s \/8li). Because we are only concerned 
with scaling, and not with the exact solution of ()20j) . a fairly crude analysis will be 
sufficient. 

First we consider the solution of (|2U|) with the noise term omitted, and the initial 
condition [see (JTHJ)] 

x (0) = (-0.7290...) (8li/cl 2 ) 1/3 . 

We define a characteristic point of the orbit, x n i(n n {), where ax 2 { ~ n n \5LL. For n < n nh 
nS/u, < dx/dn < 2n5fi, and we can approximate the noiseless orbit as 

x(n) « x(0) + a(n)(n 2 5Li), (21) 

where a(n) is a slowly varying function of n of order 1 (1/2 < a(n) < 1 for n < n n i). 
Setting ax 2 ~ n 5fi, we find that n n \ is given by 

n nl ~ (a5Li)- 1/3 , (22) 

corresponding to [c.f., Eq. (|2T]I] 

x n i ~ (<Wa) 1/3 . 

For n > n n i (i.e., x(n) > £ n i), Eq. (J2TJ|) can be approximated as dx/dn w ax 2 . Starting 
at £(n) ~ x n i, integration of this equation leads to explosive growth of x to infinity in a 
time of order (a<5/i) -1 / 3 , which is of the same order as n n \. Thus, the relevant time scale 
is (ao//) -1 / 3 [this agrees with Eq. (fTHjl in Sec. IIIEj . 

Now consider the action of noise. For n < n a \, we neglect the nonlinear term ax 2 , 
so that (|20j) becomes dx/dn = nbii + Ae(n). The solution of this equation is the linear 
superposition of the solutions of dx a /dn = n5fi and dxb/dn = Ae(n), or x{n) = x a {n) + 
Xb{n); x a {ri) is given by x a {n) = x(0) + n 2 5/i/2, and x b(n) is a random walk. Thus, for 
n < n n \, there is diffusive spreading of the probability density of x, 

AdifrH = \J (x 2 b (n)) ~ n 1/2 A. (23) 

This diffusive spreading can blur out the structure in Fig. 6. How large does the noise 
amplitude A have to be to do this? We can estimate A by noting that the periodic 
structure in Figs. 6(b,c) results from orbits that take different integer times to reach 
x ~ x n \. Thus, for n rs n n i we define a scale A n i in x corresponding to the periodicity in 
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l/6/i by [c.f., Eq. (j2U)] 



x n i ± A„i « x(0) + (n n! ± l) 2 5/z 



which yields 



A n i ~ n n i5fi. 



(24) 



If by the time n m n n \, the diffusive spread of the probability density of x becomes as 
large as A n i, then the noise starts to wash out the periodic variations with 1/Sfi. Setting 



Thus, we expect a collapse of the two parameter (A, <5/i) data in Fig. 14(a) by means of 
a rescaling of A by 5fx raised to an exponent 5/6 [i.e., A/(5fi) 5 ^ 6 ]. 

III. SCALING OF INDETERMINATE SADDLE-NODE BIFURCATIONS FOR 
A PERIODICALLY FORCED SECOND ORDER ORDINARY DIFFERENTIAL 
EQUATION 

In this section we demonstrate the scaling properties of sweeping through an indeter- 
minate saddle-node bifurcation in the case of the periodically forced Duffing oscillator 



The unforced Duffing system (i.e., fi = 0) is an example of an oscillator in a double well 
potential. It has two coexisting fixed point attractors corresponding to the two minima 
of the potential energy. For small fi, the forced Duffing oscillator has two attracting 
periodic orbits with the period of the forcing (i.e., 2tt), one in each well of the potential. 
At /I — fi* ~ 0.2446, a new attracting periodic orbit of period 6ti arises through a 
saddle-node bifurcation. In Ref. j2j]|, it is argued numerically that for a certain range 
of /i > /i* the basin of attraction of the 6n periodic orbit and the basins of attraction 
of the 27r periodic orbits have the Wada property. Thus, as /i decreases through the 
critical value the period 67r attractor is destroyed via a saddle- node bifurcation on 
the fractal boundary of the basins of the other two attractors. This is an example 
of an indeterminate saddle-node bifurcation of the Duffing system which we study by 
considering the two-dimensional map in the (x, x) plane resulting from a Poincare section 
at constant phase of the forcing signal. We consider orbits starting in the vicinity of the 
period three fixed point attractor, and, as we integrate the Duffing system, we decrease 
\i from ix s > fi* to /i/ < /i* at a small rate of 8fi per one period of the forcing signal. As 
fi approaches (with fi > fi*,) the period three fixed point attractor of the unswept 
Duffing system approaches its basin boundary, and the slowly swept orbit closely follows 
its location. For fi—fi* < small, the orbit will approximately follow the one-dimensional 
unstable manifold of the \i = /i* period three saddle-node pair. Thus, we can describe the 
sweeping through the indeterminate bifurcation of the Duffing oscillator by the theory 




0, 



'x — 0.15 x — x + x 3 = /i cost. 



(26) 
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we developed for one dimensional discrete maps. Figure 15 shows the final destination 
graph of a swept orbit initially situated in the vicinity of the period three fixed point 
of the Poincare map. The final attracting state is represented as a 1 if situated in the 
potential well where x > 0, and is represented as a if situated in the potential well where 
x < 0. As expected, the structure in Fig. 15 appears asymptotically periodic if graphed 
versus 1/Sfi. In addition to slowly sweeping the Duffing system, consider an additive 
noise term Ae(t) in the right hand side of (|26j) . where on every time step e(t) is chosen 
randomly in [—1,1], and the time step used is At = 27r/500. Figure 16(a) shows the 
dependence of the probability of approaching the attractor represented as a 1 versus the 
noise amplitude A for three specially selected values of 5fi (centers of white bands in the 
structure of Fig. 15 where the swept orbit reaches the attracting state represented by 1) 
spread over one decade. Figure 16(b) shows collapse of the data in Fig. 16(a) to a single 
curve when the noise amplitude A is rescaled by (<5/i) 5//6 , as predicted by our previous 
one-dimensional analysis (Sec. I11F|) . Thus, we believe that the scaling properties of the 
indeterminate saddle-node bifurcation we found in one-dimensional discrete maps are 
also shared by higher dimensional flows. 

IV. INDETERMINACY IN HOW AN ATTRACTOR IS APPROACHED 

In this section we consider the case of a one dimensional map f ^ having two attractors 
A and B, one of which (i.e., A) exists for all fi G \jJL s ,fjLf\. The other (i.e., B) is a node 
which is destroyed by a saddle-node bifurcation on the boundary between the basins of 
A and B, as ji increases through //* (/i* G [/i s , ///]). When an orbit is initially on B, 
and \i is slowly increased through /i*, the orbit will always go to A (which is the only 
attractor for /i > /i*). However, it is possible to distinguish between two (or more) 
different ways of approaching A. [In particular, we are interested in ways of approach 
that can be distinguished in a coordinate-free (i.e., invariant) manner.] As we show in 
this section, the way in which A is approached can be indeterminate. In this case, the 
indeterminacy is connected with the existence of an invariant nonattracting Cantor set 
embedded in the basin of A for fi > /x*. 

As an illustration, we construct the following model 

f^x) = -n + x-3x 2 -x 4 + 3.6x 6 -x 8 . (27) 

Calculations show that / M satisfies all the requirements of the saddle-node bifurcation 
theorem for undergoing a backward saddle-node bifurcation at x* = and /x* = 0. 
Figure 17(a) shows the graph of versus x at /i = /J*. Figure 17(b) shows how the 
basin structure of the map / M varies with the parameter [i. For positive values of 
fn has only one attractor which is at minus infinity. The basin of this attractor is the 
whole real axis. As /i decreases through /i* = 0, a new fixed point attractor is created at 
x* = 0. The basin of attraction of this fixed point has infinitely many disjoint intervals 
displaying fractal features [indicated in black in Fig. 17(b)]. This is similar to the blue 
basin B\p] of the attractor of the previous one-dimensional model (see Sec. Ill A]) . 

The blue region in Fig. 18(a) is the basin of attraction of the stable fixed point 
destroyed as \i increases through //*. For every value of /i we consider, the map has 
invariant Cantor sets. The trajectories of points which are located on an invariant Cantor 
set, do not diverge to infinity. One way to display such Cantor sets, is to select uniquely 
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defined intervals whose end points are on the Cantor set. For example, Fig. 18(a) shows 
green and red regions. For every fixed parameter value p, the collection of points that 
are boundary points of the red and green regions, constitutes an invariant Cantor set. In 
order to describe these green and red regions, we introduce the following notations. For 
each parameter value p, let p M be the leftmost fixed point of / M [see Fig. 17(a)]. For every 

xo < Pn, the sequence of iterates {x n = /jf^o)} is decreasing and diverges to minus 
infinity. For each value of p, let be the fixed point of / M to the right of x = at which 
7fcf(9/J > 1- A point (x; p) is colored green if its trajectory diverges to minus infinity 
and it passes through the interval (g M , oo), and it is colored red if its trajectory diverges 
to minus infinity and it does not pass through the interval (g M , oo). Denote the collection 
of points (x\ p) that are colored green by G[p], and the collection of points (x;p) that 
are colored red by R[p]. Using the methods and techniques of [i^, it can be shown that 
the collection of points (x\ p) which are common boundary points of G[p] and R[p] is 
a Cantor set C[p] 0]. In particular, the results of [22| imply that for p = //* = the 
point x* = belongs to the invariant Cantor set Cf/i*]. 

Figure 18(b) is a zoom of Figure 18(a) in the region of the saddle-node bifurcation. 
For values of p > /i*, in the vicinity of (x*; one notices a fractal alternation of red 
and green stripes. The green and red stripe structure in Fig. 18(b) shares qualitative 
properties with the structure in Fig. 2(b). All the analysis in Sec. HU can be adapted 
straightfowardly to fit this situation. 

Figure 19 shows how the chaotic saddle of the map / M varies with p. The chaotic saddle 
is generated numerically using the PIM-triple method. For an explanation of this method 
see Nusse and Yorke 0]- Using arguments similar to those in Sec. Ill C[ we predict that 
changing the horizontal axis of Fig. 19 from p to (p — p*)^ 1 / 2 makes the chaotic saddle 
asymptotically periodic. Numerical results confirming this are presented in Fig. 20(a). 
For fp given by (|27|). we were able to find a parameter value p** = 0.23495384 where 
changing the horizontal axis of Fig. 19 from p to (p — /i**) -1 / 2 [see Fig. 20(b)] apparently 
makes the chaotic saddle asymptotically periodic [with a different period than that of 
Fig. 20(a)]. As in the case discussed in Sec. HU past the saddle-node bifurcation of / M 
at p*, infinitely many other saddle- node bifurcations of periodic orbits take place on the 
invariant Cantor set C[p]. We believe that p** is an approximate value of p where such 
a saddle-node of a periodic orbit takes place. 

V. DISCUSSION AND CONCLUSIONS 

In this paper, we have investigated scaling properties of saddle-node bifurcations that 
occur on fractal basin boundaries. Such situations are known to be indeterminate in 
the sense that it is difficult to predict the eventual fate of an orbit that tracks the pre- 
bifurcation node attractor as the system parameter is swept through the bifurcation. 
We have first analyzed the case of one- dimensional discrete maps. Using the normal 
form of the saddle-node bifurcation and general properties of fractal basin boundaries, 
we established the following universal (i.e., model independent) scaling results 

• scaling of the fractal basin boundary of the static (i.e., unswept) system near the 
saddle-node bifurcation, 

• the scaling dependence of the orbit's final destination with the inverse of the sweep- 
ing rate, 
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• the dependence of the time it takes for an attractor to capture a swept orbit with 
the -1/3 power of the sweeping rate, 

• scaling of the effect of noise on the final attractor capture probability with the 5/6 
power of the sweeping rate. 

All these results were demonstrated numerically for a one-dimensional map example. 
Following our one-dimensional investigations, we have explained and demonstrated nu- 
merically that these new results also apply to two-dimensional maps. Our numerical 
example was a two-dimensional map that results from a Poincare section of the forced 
Duffing oscillator. In the last section of the paper, we have discussed how the new 
results listed above apply to the case where a saddle-node bifurcation occurs on an in- 
variant Cantor set which is embedded in a basin of attraction, and we have supported 
our discussion by numerics. 

This work was supported by ONR (Physics), NSF (PHYS 0098632), NSF (0104087), 
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FIG. 1: Construction of the function fu{%) starting with (a) the third iterate of the logistic map, 
g{x) = rx{l — x), with r = 3.832, and adding a perturbation (b) /isin(37rx) (p = 5.4 x 10~ 3 ). 



FIG. 2: (a) Basin structure of the map versus the parameter \x on the horizontal axis 
(0 < fi < 5.4 x 10~ 3 and < x < 1). The attractor having the blue basin is destroyed 
at (j, « 2.79 x 10~ 3 . (b) Detail of the region shown as the white rectangle in Fig. 2(a), 
2.75 x KT 3 < n < 3.55 x 1(T 3 and 0.145 < x < 0.163. 



FIG. 4: (a) Detail of Figure 2(b), with the horizontal axis changed from fi to (fi — fi*)^ 1 ^ 2 
for /i > /x*; 2.75 x 10~ 3 < fi < 3.55 x 10~ 3 and 0.145 < x < 0.163 The green stripes from 
Figure 2(b) are colored black and the red stripes are colored white. The approximate position 
of the point x* where the saddle-node bifurcation takes place is shown. x c indicates the nearest 
critical point, (b) Detail of Figure 3, displaying how the box dimension D of the fractal basin 
boundary varies with l/(^* — /u) 1 / 2 . The horizontal axis of Figures 4(a) and 4(b) are identical. 



FIG. 5: Qualitative graphs of the solution of Eq. (JSJ), jl n (xo), for three consecutive values 
of n. Note the horizontal asymptotes [/i -1 / 2 = (n — l)a 1 / 2 7r, na 1 / 2 ^, and (n + l)a 1 / 2 7r], the 
vertical asymptotes [x s = (a(n — 1)) , (an) -1 , and (a(n + 1)) ] 3 both shown as dashed lines, 
and the intersections of the solid curves with xq = which are marked with black dots. 

FIG. 6: (a) Final attracting state of swept orbits versus We have chosen fi s = ji s + fi* = 0, 
and fif = 4.5 x 10~ 3 . The attractor R flf is represented by 1 and the attractor G^ f is represented 
by 0. (b) Detail of Fig. 6(b) with the horizontal scale changed from Sfi to l/5ji. The structure 
of white and black bands becomes asymptotically periodic, (c) Final state of orbits for the 
system fa versus 1/5 fi. The final state of an orbit is defined to be if there exists n such 
that 100 < x n < 250, and is defined to be 1, otherwise. We have chosen jl s = — fi*, so that 
Figs. 6(b,c) have the same asymptotic periodicity. 



FIG. 3: Fractal dimension of the basin boundary versus [x. Notice the continuous variation 
for n < fi* and the discontinuous jump at fi*, the parameter value at which the saddle- node 
bifurcation on the fractal basin boundary takes place. 
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FIG. 7: Numerical results for the inverse of the limit period in 1/Sfj, versus fi s . The fit line is 
[A (1 / 6 /J-)] -1 = —0.9986^ s +0.0028 and indicates good agreement with the theoretical expla- 
nation presented in text. 



FIG. 8: Graphs of fp,(x) at different values of the parameter fx. The black dots indicate the 
stable fixed points of for different values of fi. 



FIG. 9: The calculated fractal dimension D' of the structure in the intervals between the centers 
of consecutive wide white bands in Fig. 6(b) versus their center value of l/8ji. 



FIG. 10: Capture time by the fixed point attractor G iif versus 1/Sfj,. We have chosen ji s = 0. 
The range of is approximatelly one period of the graph in Fig. 6(b), with 5fi ~ 10~ 8 . The 
vertical axis ranges between 250 and 650. No points are plotted for values of 5fi for which the 
orbit reaches the fixed point attractor R fMf . 



FIG. 11: Capture time by the middle fixed point attractor of / M versus 5/j, {fx s = 0). The best 
fitting line (not shown) has slope -0.31, in agreement with the theory. 



FIG. 12: Probability that one orbit reaches the middle fixed point attractor of versus the 
noise amplitude A, for five different values of 6/j, (10^ 5 , 10~ 5 ± 2.5 x 10~ 8 and 10~ 5 ± 5 x 10~ 8 ). 
We have chosen fi s = 0. 



FIG. 13: Probability that an orbit reaches the middle fixed point attractor of f^, for five 
selected values of 5[i spread over two decades: (a) versus the noise amplitude A, and (b) versus 
A/(5[i) 5 / 6 , We have chosen fi s = 0. From right to left, the S/j, values corresponding to the 
curves are approximately: 3.445974xl0" 5 , 1.147767xl0" 5 , 3.820744x 10~ 6 , 1.273160x 10~ 6 
and 4.243522xl0~ 7 . 
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FIG. 14: Probability that an orbit of fp, reaches a fixed interval far from the saddle-node 
bifurcation (i.e., [100, 250]), for five values of Sfi spread over two decades: (a) versus the 
noise amplitude A, and (b) versus A/(5[i) 5 / e . We have chosen fi s = 0. From right to left, 
the 5/i values corresponding to the curves are approximately: 3.451540xl0~ 5 , 1.149162xl0~ 5 , 
3.829769xl0~ 6 , 1.276061xl0~ 6 and 4.253018xl0~ 7 . 



FIG. 15: Final attracting state of swept orbits of the Duffing oscilator versus 1/Sfi. The 
structure of white and black bands becomes asymptotically periodic. We have chosen /j, s = 
0.253, and Hf = 0.22. The attractor in the potential well for x > is represented as a 1, and 
the attractor in the potential well for x < is represented as a 0. 



FIG. 16: Probability the Duffing oscillator reaches the attracting periodic orbit in the po- 
tential well at x > for three values of 5(jl spread over one decades: (a) versus the noise 
amplitude A, and (b) versus A/(5n) 5 ^. We have chosen fi s = 0.253. From right to left, the 
Sfi values corresponding to the curves are approximately: 4.628716x 10™ 5 , 1.461574x 10~ 5 and 
4.621737xl0~ 6 . 



FIG. 17: (a) Graph of /«(#) versus x at the bifurcation parameter, (b) Basin structure of map 
ffi{x) versus the parameter ^ (—0.3 < n < 0.3 and —2 < x < 2). The basin of attraction of the 
stable fixed point created by the saddle-node bifurcation is black while the basin of attraction 
of minus infinity is left white. 
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FIG. 18: (a) Basin structure of / M versus fi (-0.3 < [i < 0.3 and -2 < x < 2). We split the 
basin of attraction of minus infinity into two components, one plotted as the green region and 
the other plotted as the red region. The green region is the collection of all points that go to 
minus infinity and have at least one iterate bigger that the unstable fixed point q^. The red 
set is the region of all the other points that go to minus infinity, (b) Detail of Fig. 15(a) in the 
region shown as the white rectangle, —0.005 < ji < 0.015 and —0.09 < x < 0.41. 
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FIG. 19: The chaotic saddle of / M versus \i (—0.3 < \i < 0.3 and — 2 < x < 2) generated by the 
PIM-triple method. 



FIG. 20: The chaotic saddle of the map in the vicinity of the saddle-node bifurcation with 
the horizontal axis rescaled from fi to: (a) — fi)' 1 / 2 . Notice that the chaotic saddle becomes 
asymptotically periodic (-0.008 < x < 0.337, 10 < (/x*-/x) -1 / 2 < 15). (b) (ii^-fi)' 1 / 2 , where 
= 0.23495384. We believe that corresponds to the approximate value of the parameter 
H where a saddle-node bifurcation of a periodic orbit of / M takes place on the Cantor set C\p]. In 
this case, the chaotic saddle also becomes asymptotically periodic (—0.162 < x < 0.168, 9.97 < 
(M** - /u)- 1 / 2 < 2010). 
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